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Abstract 

We show that for any semi-random transposition shuffle on n cards, the mixing time 
of any given k cards is at most nlogfc, provided k = o((n/ log n)^/^). In the case of 
the top-to-random transposition shuffle we show that there is cutoff at this time with a 
window of size 0(n), provided further that fc ^ oo as n ^ oo (and no cutoff otherwise). 
For the random-to-random transposition shuffle we show cutoff at time (l/2)rtlogA: for 
the same conditions on k. Finally, we analyse the cyclic-to-random transposition shuffle 
and show partial mixing occurs at time < an log A: for some a just larger than 1/2. We 
prove these results by relating the mixing time of k cards to the mixing of one card. Our 
results rely heavily on coupling arguments to bound the total variation distance. 

1 Introduction 

We denote by Sn the symmetric group on n elements which we shah view as a deck of n 
labeled cards. We use [n] to denote {1, . . . ,n}. We are interested in a notion of mixing of 
these n cards when the shuffling mechanism is semi-random transposition shuffling. At every 
step of such a shuffle we choose 2 cards - one with our left hand and one with our right 
(independently of each other and of the past choices). We then switch the positions of the 
two selected cards. We denote by Lf and Rt the location chosen by the left hand and right 
hand at time t, respectively. The left hand is allowed to choose a card according to any rule, 
deterministic or stochastic, and it may also depend on time. The right hand chooses a card 
uniformly. 

As well as stating a result for general semi-random transposition shuffles, we also focus on 3 
particular shuffles. These are the cyclic-to-random transposition shuffle, with Lt = t mod n, 
the top-to-random transposition shuffle, with Lt = 1, and the random-to-random transposi- 
tion shuffle, with Lt ~Unif([n]). 

We first state our results and then present some background and motivation. We shall denote 
by Aici,...,Cfc the joint law of the locations of cards ci, . . . ,Cfc in a permutation a. Let at be 
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a permutation at time t which starts at time from permutation do and evolves by a given 
semi-random transposition shuffle. Let vr be a uniform permutation. We define dci,...,Cfe(i) by 



4i,...,cfe(i) := max \\^lcl,...,Ck - /^ci,...,c^ 



TV, 



where ||// — z^||tv is the total variation distance between two measures fi and v. Let denote 
the set of all subsets of [n] of size k. 

Definition 1. We define the k-partial mixing time (at level e) of a shuffle to be 
4ix(e) := minjt > : max dc,...,cdt) < e|. 



Our main result is the following. 

Theorem 1. Suppose k = o{{n/ log n)^^^) and fix e > 0. For any semi-random transposition 
shuffle with k-partial mixing time tmix(^) '^'^'^ ^''^V ^ > there exists no = no{6) such that for 
all n > uq, 

{i). imix(^) ^ ^mix((£ ~ ^)/^)' '^^'^ moreover, 
(n). 4,(1/4) <n(log A: + 3/2). 



Theorem [T] allows us to calculate an upper bound on the /c-partial mixing time of any semi- 
random transposition shuffle by just considering the movement of one card in the deck. 

Theorem 2. Suppose k = o((n/ log n)^/^). For the cyclic-to-random transposition shuffle, 
there exists a constant C > 0, such that for all n sufficiently large, 

4i,(l/4) < 0.5006n(log A: + C). 

Definition 2. We say that a shuffle has k-partial cutoff at time tmix with a window of size 
oj{n) = o(tmix) if the following two conditions hold: 

(i). lim limsup max dci Cfe(imix + a'^('^)) = 0, 

a^oo „^oo (ci,...,Cfc)Gnfc 

[a). lim liminf max dc^ c (tmix + aa;(n)) = 1. 

a-s--oo n^oo (ci,...,Cfc)GQfe 

Theorem 3. Suppose k = o((n/ log n)^/^) and also A; — )• oo as n ^ oo. Then the top-to- 
random transposition shuffle has k-partial cutoff at time nlogk with a window of size 0{n). 
Furthermore, if k < K , for some constant K for all n, then there is no cutoff. 

Theorem 4. Suppose k = o((n/ log n)^/^) and also k —?■ oo as n ^ oo. Then the random- 
to-random transposition shuffle has k-partial cutoff at time 0.5nlog/c with a window of size 
0{n). Furthermore, if k < K, for some constant K for all n, then there is no cutoff. 



The study of mixing times of Markov chains and the search for cutoff is a much-studied area 
of probability. In terms of walks on the symmetric group, random transpositions are one of 
the most natural and simplest models. It was first shown bv lDiaconis and Shahshahanil 1981 ] 
that the random-to-random transposition shuffle has cutoff at time O.Snlogn with a window 
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of size 0{n). Their technique uses Fourier analysis on the symme tric group. Since then, this 
result has been shown with a strong-sta, t ionar y time argument by Matthews _ 1988ll an d with 
couphng arguments by Berestvcki et al. 2011 ] and independently by BlumberX 2011 ] . 



Theorem [T] gives an upper bound of n log k on the /c-partial mixing time of any semi-random 
transp osition shuffle. Regard ing full mix i ng, it has been shown by ISaloff-Coste and Zuiiiga 



2007l | and independently by iGanapathvl 2007l | that the mixing time of any semi-random 



transposition shuffle is at most nlogn, with the top-to-random having cutoff at this time. 



The cyclic-to-random transposi tion shuffle was invented by iThorpI 1965l | and the question 
of its mixing time was posed bv lAldous and DiaconisI 1986t |. Unlike with the top-to-random 
and random-to-random shuffles, after just n steps of the cyclic-to-random shuffle, every card 
has almost surely been selected with either the left hand or the right hand at least once. 
The standard coupon-collector argument for establishing a lower bound does not give the 
correct order for the mixin g time of this shuff le. The current best known lower bound is 
about 0.12nlogn, shown bv lMossel et al.l 2004l | obtained by analysing the eigenfunctions of 
the transition matrix for the movement of a single card. The best known upper bound is 
the general n log n of any semi-random transposition shuffle. It remains an open question to 
determine if and when cutoff occurs for this shuffle. 

The notion of studying the evolu tion of only some of the cards in a deck has been previously 



considered by lAssaf et al.l 20111 ]. Here they calculate the distribution of the location of 
one card after a step of riffle shuffling. For the cyclic-to-random transposition shu ffle the 



distri bution of the location of one card aft er a round (n steps) has been calculated by iPinskv 



2012l | and independently bv lPvmarl 201l|. It is a natural extension to try to understand how 
k cards in a deck evolve and our results reveal interesting dynamics. As with full mixing, 
the top-to-random shuffle is the slowest shuffle in terms of partial mixing, and using the 
random-to-random shuffle halves this time. We also find that the cyclic-to-random does not 
take much longer to mix k cards than the random-to-random (and may even be much faster). 

Structure of the rest of the paper: In section[21 we give the proof of Theorem[T]by first studying 
the movement of one card and then the joint movement of a given k cards. Next, in section [3] 
we show partial cutoff of the top-to-random and random-to-random shuffles by using Theorem 
dJ We then give our upper bound of the /c-partial mixing time of the cyclic-to-random shuffle 
in section [H before finally presenting some open questions in section [5j 



2 Proof of Theorem [T] 

Definition 3. For each non-negative integer t, let it be a distribution on [n]. We say that a 
permutation evolves by {lt)t>{)-to-random if it evolves according to a semi-random transpo- 
sition shuffle and for each t the choice of the left hand at time t is chosen according to the 
distribution If 

We begin by studying the mixing time of one card for any (^t)j>o-to-random transposition 
shuffle. Let (o"t)t>o denote the state at time t of a permutation that evolves by (^t)t>o-to- 
random, started from permutation a^, and (vrj)t>o a permutation also evolving by {it)t>o-^o- 
random, started from a uniform permutation ttq. We denote by o""^(i) the location occupied 
by the card with label i in a permutation a. Recall that for each t, Lt and Rt are independent 
from each other and from past choices of the left and right hands. 
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Lemma 1. For every {£t)t>o-to-random transposition shuffle, 

maxdi(t) := max max 11//°^* ~y^7*||TV < exp(— t/n). 

ie[n] ie[n] <^o 

Proof. We use the coupling description of the total variation distance, that is, 
\\jjL — z^IItv = inf{P(X ^ Y) : {X, Y) is a coupling of jj, and v}. 

For each we use the same Lt for both processes at and -Kt- We denote by t— the time just 
before the t^^ transposition. 

To describe the joint evolution of at and ttj, we choose at time t the card in location Rt 
with the right hand in both processes unless Rt € {a'[_}{i),'K^}{i)}. If Rt = a'[_}{i), we select 
position a'[_}{i) in deck vr^ and select position 7r^^(i) in deck at (so that card i is not chosen in 
either deck). On the other hand, if Rt = T^t-i^) select position i^t-i^) and position 

in at (so that card i is chosen in both decks). It is clear that the marginal distributions 
are correct - in each deck the probability of selecting a particular card with the right hand is 
1/n, independently of all past choices. 

The result of this coupling is that once card i is chosen with the right hand in (vr^), the 
locations of card i in both processes will become and then stay forever equal. The time to 
select card i in (ttj) is Geom(l/n). Thus we have 

P((7,-^(i) = ni'ii)) > 1 - expi-t/n), 

and the bound on the total variation distance follows. □ 

We refer to the k cards of interest ci, . . . ,Cfc as the special cards. We present now a bound 
on the number of times the left hand chooses one of the special k cards up to time t, for 
any (^t)t>o-to-random transposition shuffle. Choose such a shuffle and let {Lt)t>o be an 
independent sequence of random variables such that Lt ~ it for each t. 

Lemma 2. There exists a constant C such that for any choice of k special cards, for each 
t > and n sufficiently large, 

s=l 

Proof. First note that we can write 

t t k k t 

s=l s=l i=l 1=1 s=l 

Let Ni{t) be the number of times that the left hand selects card i by time t, that is, 

t 

s=l 
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Let Ti{m) be the time between the (m — 1)**^ and the m}^ selection of card i with the left 
hand, that is, 

m—1 m—1 

Ti{m) := m.m{t: t>Y^ Ti{£), Lt = (Tt_}ii)} - Ti{i). 

e=i t=i 

Then Ni{t) satisfies 

X 

{Niit)>x} = {Y,n{m)<t}. (1) 

m=l 

Notice now that Ti{m) > fi{m) where fi{m) is the amount of time after the (m— 1)*^ selection 
of card i with the left hand and the next time card i is selected (with either the left or right 
hand), that is, 

m—1 m—1 

n{m) := mm{t:t>Y, ^^W' ^t-'W e {Lt,Rt}} - ^ Ti{£). 

i=i e=i 

We further define Tj(r) to be the time at which card i is chosen with the left hand for the r*^ 
time. That is, Ti{r) = Ylm=i'''ii''^)- Note that for s satisfying 

T,{r) <s <T,ir) + n{r + l), 

wc have cr'^}{i) = RT.-{r)i where Rt is the location chosen with the right hand at time t. We 
show that for each m > 1, E(fi(m)) > Cin, for some constant ci. 

Firstly, at each time s > Ti{m — 1) we have probability 1/n of choosing card i with the 
right hand. On the other hand, to select card i with the left hand it has to be in a location 
accessible by the left hand. For example, in the case of top-to-random, the left hand can only 
ever select the card at the top of the deck. Therefore in this case if at time t card i was in 
any position other than at the top of the deck, it would certainly be chosen with the right 
hand before the left hand. Furthermore, since after the left hand selects card i it is moved to 
a uniform location, in order to minimize the expected time to select this card again, the next 
n positions selected by the left hand Lrp.(^m-i)+ij ■ ■ ■ )-^Ti(m-i)+n should be a bijection of [n]. 
We formalize this idea below. We let 

H = mi{h > 1 : LT.(^jn-i)+h = -^ri(m-i)}) 

i? = inf{r > 1 : RTi(m-l)+r = ^Ti(m-l)}- 

Note that H and R are independent and that R is geometrically distributed with success 
probability 1/n. Further, in the case where Lj^.(^jn_i^^i, . . . , Lrp.^^_-^'^_^^ is a bijection of [n], 
H is uniformly distributed on [n]. 

Depending on the semi-random transposition shuffle we axe using, the {fi{m))m>i may not 
be independent. However, we claim that for each m > 1, -fj(m) stochastically dominates 

fi{m), independent identically distributed random variables with distribution mm{H, R). We 
use a coupling argument to show this stochastic domination. We define fi{m) to be the 
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amount of time after Ti{m — 1) until the left hand or right hand chooses card i but now 
with different L- values. We denote these new n L- values after time Tj(m — 1) as Lj{m — 1), 
^^{m — 1), . . . ,L2{m — 1). This list is a bijection of [n] with the property that if the first 
occurrence of the value j in L-^.i-^.i).,.!, Ly.(^_i)_|_2 • • • , appears before the first occurrence of 
j', that is if 

min{i > : LT.(^_i)+t = j} < mm{t > : = j'}, 

then j must appear before / in Lj{m — 1), Lf{m — 1). 

We let J = {j G [n] : 3t > : Lrp.(j^_i^^t = j}. For j G [n\ which do not occur in 
-LTi(m-i)+i) -^Ti(m-i)+2) • • •) j G N \ J ■, we use the rule that they are listed arbitrarily in 
L^(m — 1), . . . , L"(m — 1) but after all values that do occur. That is, 

Vji G J,i2 e N \ J, if L*i = ^ thenti < is. 



At time Tj(m— 1) card i is moved to location R^^i^^-i) and will remain there until chosen with 
the left hand or the right hand. The choice of Lj{m — 1) for j > 1 ensures that fi{m) < fi{m). 
The independence and identical distribution properties of fi{m) for m > 1, come from the 
fact that the location card i is put into every time it is chosen with the left hand is chosen 
independently and uniformly on [n] and that L|(m — 1), . . . , L"(m — 1) is a bijection of [n]. 
Taking fj(l) = 0, we have that 

X X 

m=l m=l 



It follows that for m > 1, 

E(fi(m)) = -E(fi(m) \ H = h) = -J^Yl \H = h,R = r)V{R 

h=l ^ ^ h=l r=l 

^ n n 

= -J2Yl ^(^^M \H = h,R = r){l- 1/ny-' - 

h=l r=l 

n h—1 n 

= E ( E - (1 - vn)^^^ + E Ml - i/nr' ) 



h=l r=l r=h 



^(n-3)(l-l/n)" + l>fein, (3) 



for some constant ki. 
Similarly, for each m > 1, 

n h—l 



/i=l r=l r=h 
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for some constant k2- Using this and equation ([3]), we have that Var('rj(m)) = k^n'^, for some 
constant k^. Since the Tj(m) are independent for all m > 1, we deduce that 

X 

Var(^ fi{m)) = ksix - l)n^ . (4) 

m=l 

From equations ([T]) and ([2]), we have that 

t k k t 



s=l 1=1 1=1 x=l 



i=l x=l m=l 1=1 x=l m=l 

By Chebyshev's inequality and equation dU, for any a > 0, 

m=l m=l 

Using equation ([3]) this gives us 



m=l 



and thus for t <ki{x — l)n, 



TO=1 

We now have 



{ki{x - 1) - t/n) 



2 ■ 



t X \l+t/(kin)'] X t X 

x=l m=l x=l m=l x=^i^t/(^kin)]+l m=l 

t X 

<2 + t/{kin)+ Y IP(j^^*M<i) 

a;=[l+t/(fcin)]+l m=l 

<2 + t/(A;in)+ V -— ^^^^ —-^r 

a;=ri+t/(fcin)]+l ^ ^ ' ^ 

< 2 + t/(A;in) + k^ logi 

< c(- + logt 

for some constants /c4,C. Using this and equation 1^ completes the proof. □ 

Let it|, . . . , o"^ denote the state at time t oi k permutations, each evolving independently and 
by (^f)t>o-to-random shuffle and starting from a permutation ctq- We describe a coupling of 
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Figure 1: The desired property of the coupling 



these k permutations and a (A;+ 1)*^ permutation denoted which starts from permutation 
(To and also evolves by (^t)t>o-to-random shuffle (but not independently from the other per- 
mutations). The coupling has the property that initially and for a period of time afterwards, 
the locations of special card Cj in permutation o"* for each \ <i <k matches the location of 
card Cj in permutation ^ see Figure [TJ 

Lemma 3. There exists a constant c such that for any choice of k special cards, for each 
t > and n sufficiently large, 



I 



fJ'C 



X /^cfe IItv 



< c 



tk^ fc^ log t 
1 — 



n 



Proof We use a coupling argument to bound this total variation distance. We use the same 
choices of Li, L2, ■ ■ for all permutations. For the right hand, we let Rl denote the location 
chosen by the right hand at time t of permutation a^, 1 < i < k (which are independent and 
uniformly distributed on [n]). We describe how to choose R^, the choices of the right hand in 
permutation a^, so that it evolves by (^t)t>o-to-random with the desired property of Figure 
[1] for as long as possible. 

Firstly, if Lt selects one of the special cards in permutation a^, say Lt = ((T^_)~^(cj), then we 
set R^ = R\. On the other hand, if Lt does not select one of the special k cards, then we toss 
a coin which lands heads with probability p, independently on each toss, where 



p :-- 



k/n + (1 - l/nY 



> 1 - c/cV"^ 



(6) 



for some constant c. If the coin lands tails we choose U uniformly on {((Tj_)~^(ci), 1 < i < k} 
and set R^ = U, that is, we choose a special card, each with equal probability. If the coin 
lands heads, we choose U uniformly on [n]. If U chooses a special card, that isU = ((Tf_)~^(cj) 
for some 1 < i < /c, we set i?? = Rl IfU chooses a non-special card we set R^ = U unless at 
least one of the permutations cr^, . . . , ct'^ selects its special card with the right hand, that is 



E 



{l<i<k: R\ = {al_)-'{c,)}^^ 
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In this case we choose i uniformly from the set E and let = R\ = (cr|_) ^(q). 

We show that at every step t the choice of the right hand i?^ is uniformly distributed on 
[n]. Firstly, conditioned on the event {Lt = {a^_)~^{ci), some 1 < « < /c}, then the choice of 

is clearly uniform. We have to check more carefully for the case where we condition on 
the complement of this event. We fix non-special card j, and calculate the probability that 

= We have 

= ^(^ + (l-l/n)M=l, 
n\n / n 

as desired. It follows by symmetry that the probability of selecting with the right hand a 
specific special card is also 1/n at every step. 

We now let At be the event that up to time t, the location of card Cj in permutation cr* is 
equal to the location of card Cj in permutation a^, for all 1 < i < k, that is, 

At := = (a°)-^(ci), for aU < s < t, 1 < i < k}. 

Certainly Aq holds. We list below the various situations which can result in a mismatch of 
the locations of the special cards at time t. If none of these situations occur by time t then 
event At will hold. 

1. Suppose Lt = (o"['_)^^(cj), for some 1 < i < k. A mismatch if Rl = (o"|_)^^(cj') for 
some i' ^ i, which happens with probability (k — l)/n, or i?f = (o"j!_)^^(cj') for some 
i' 7^ i, which happens with probability 1 — (1 — l/n)'^~^. 

2. Suppose Lt chooses a non-special card, the coin lands tails and we choose some permu- 
tation, say fj*, to copy for the right hand. A mismatch if Rl = {al_)^^{ci') for some 
i' 7^ i, which happens with probability {k — l)/n, or Rl = {al_)~^ {ci') for some i' ^ i, 
which happens with probability 1 — (1 — l/n)^~^. 

3. Suppose Lt chooses a non-special card, the coin lands heads and U chooses a special 
card, say Cj. A mismatch if Rl = ((T^_)~^(cj') for some i' ^ i, which happens with prob- 
ability {k — l)/n, or Rt = (o"j!_)^^(cj/) for some i' ^ i, which happens with probability 
1 - (1 - l/n)^-^. 

4. Suppose Lt chooses a non-special card, the coin lands heads and U chooses a non- 
special card. A mismatch if |{1 < i < A; : Rl = (<7j_)^^(cj)}| > 1, which happens with 
probability 1 — (1 — 1/n)'^ — {k/n){l — 

For 1 < J < 4, we let Ej{t) denote the number of times situation j above occurs by time t. 
We use Lemma [2] to control Ei{t). Let C be the constant in Lemma [2] and c the constant in 
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equation ([6]). We have 

1.2 

E(Ei(t)) < Ckit/n + logt)((k - l)/n + 1 - (1 - l/n)^-^) < 2C—(t/n + logt), 

n 



E{E2{t)) < t{l - p){{k - l)/n + 1 - (1 - l/n)"-^) < 2ct^, 



E{E3{t)) < t-{{k - l)/n + 1 - (1 - l/n)'^-!) < 2t^, 

E{Ei{t)) < t{l - (1 - 1/nf - {k/n){l - < £. 

We therefore have by Markov's inequahty, 

P(At) > P(^i(t) + E2{t) + E^{t) + E^{t) = 0) 

> 1 - nEi{t)) - n^t)) - E(^3(t)) - nEiit)) 

/tk'^ k'^logt\ 

>l-c — + ^ , 

V n / 

for some constant c. This completes the proof. 

Putting together Lemmas [1] and [3] we are able to prove Theorem [TJ 



□ 



Proof of TheoremlJl Let ttq be a uniform permutation. Let vr^, . . . , ttj^ denote the state at time 
t oi k permutations, each evolving independently by (£f)t>o-to-random shuffle and starting 
from permutation ttq. Further, let vrf be another permutation at time t evolving by the same 
shuffle and also starting from ttq (however, not evolving independently from the others). 

By the triangle inequality for total variation we have 

II ''t TTO It II <^t II / II 4 <^t\\ 

, II f^f "t '^t '^t 1 1 

+ ll/ici* X ■■■ X fj.cl - X ■■■ X ficl IItv 

+ ll/^ci* X ■■■ X fl4 - /ici*,...,cJ|TV- 

For the first and last inequalities we use Lemma [3] to give 

\\f^cl...,c, - ^^cl...,cJ\^v < 2c(^^ + — - — ) + Wf^cl X-- - xficl -f^cl X • • • X ficl IItv 

<2c(^ + —)+k max /i,'-/xr°TV- 

Maximizing over the choice of the k special cards and noticing that the first term on the 
right-hand side is o(l), we obtain for all n sufficiently large, 

max ||/icL..,Cfe - /ici,...,cJlTV < (5 + A; max 11//.* -/if°||TV- 
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This holds for any do and thus as n — > oo, 



max dci,...,Cfc(i) < 0(1) + A;max(ij(t). (7) 

(ci,...,Cfc)enfc ' ' ie[n] 



We deduce that 



4ix(£) := min{t > : max ||^S,...,Cfc - /^ci,...,cJItv < e} 

(ci,...,Cfc)Gf2fc 



< min{t >0:6 + /cmax — A^r^llTV < e} 

«G[n] 



min{t > : max ||^^* - ^f^TV < - 5)/^} 



iG n 



= ti,,{{e-6)/k). 

We apply Lemma [1] which says that i^ixC^) — — "-loge to deduce that tmix((^ ~ ^)/^) ^ 
-nlog((e - Thus imix(l/4) < n(logA; - log(l/4 - 5)) and taking 5 = 1/4- e"^/^ gives 

the desired result that t^^Jl/^) < n(log k + 3/2). □ 



Remark 1. The main idea of this proof is to show that the movement of the k special cards 
is close to independent. Our restriction on the value of k comes into play here - for k larger 
than n^/"^ their movement will in fact no longer be close to independent (in the sense that 
there will be times at which the left and the right hands choose a special card at the same 
time). 



3 Cutoff for top-to-random and random-to-random 

We begin this section by showing cutoff of the partial mixing time of the top-to-random 
transposition shuffle. The lower bound is essentially the coupon-collector problem. 

Proof of Theorem O We shall show that for the top-to-random transposition shuffle the fol- 
lowing conditions hold: 

(z). limsup max dci,...,Cfc("'log A; -|- an) < e~". 

n— >oo (ci,...,Cfc)Gnfe 

(ii). Suppose /c — )• 00 as n — )• 00. Then 

lim liminf max dc-, Cfc('T'log/c + an) = 1. 

00 ra-i-oo (ci,...,Cfe)GOfe 

[Hi). Suppose k is bounded above by constant K for all n. Then 
liminf max dci ci {n\og k + an) > e'^"^ / K. 

n->-oo (ci,...,Cfc)Gnfc 

Prom equation ([7]), we have 

max (ici,...,Cfc('^logA; + an) < o(l) -|- /c max dj (n log A; -|- an). 

(ci,...,Cfc)Gr2fc ' ' iG[n] 

However, from Lemma [H 

max dj (n log + an) < e^"/fc. 

iG[n] 
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For statement (ii), we show that for each £ > 0, 

hm hminf max dc, ci. ('t- log + an) > 1 — e. 

a-^-oo n->oo (ci,...,cfc)enfc ' ' 

Fix £ > and let C = 2/e. For a choice of C = {ci, . . . , c^}, let TZt = {(Ti_(i?i), . . . , at-{Rt)} 
be the labels of cards chosen by the right hand up to time t. We let T denote the first time 
that exactly C of the k special cards are yet to be selected with the right hand in the evolution 
of at, 

r = inf{t>0: \ntnc\ = k-c}. 

We let E(fj,) be the event that permutation n G Sn has more than C fixed points. Since 
{T >t} C E{at), it suffices to show that P(T > n log A; - an) > 1 - e/2 and F{E{Tr)) < e/2 
for n and then a sufficiently large for a uniformly chosen permutation tt. 

We denote by Xf the number of special cards that have not been moved by time t in the 
evolution of at- We show Xt is concentrated around its mean. We let Ai{t) be the event 
that card q has not been selected by time t. Then we can write Xt = Yli=i '^{Ai{t)}j so 
E{Xt) = k{l — l/nf. Furthermore, 

k k 

= k{l - 1/n)* + k{k - 1)(1 - 2/n)*. 

Thus we have 

Var(Xt) = A;(l - 1/n)* + k{k - 1)(1 - 2/n)* - k'^{l - 1/nf* 

< k{l - 1/n)* + k'^ [(1 - 2/n)* - (1 - l/n)2*] 

< A;(l - 1/n)*. 

Recall that /c ^ oo as n — )■ oo. By Chebyshev's inequality it follows that 

F{\Xt-E{Xt)\>E{Xt)-C) 

Var(Xt) 
{E{Xt) - Cf 

k{l - 1/n)* 
{k{\ - 1/n)* - C)2 

£/2, 

for all n and then a sufficiently large, with t = n log k—an. We are left to show that P(£'(7r)) < 
e/2 for all n sufficiently large. However, this follows trivially by Markov's inequality since the 
expected number of fixed points in a uniformly chosen permutation converges to 1 as n ^ oo 
and C = 2/e. 

For statement (iii), we let E{ii) be the event that special card ci is in location aQ^{ci) in 
permutation ji. Clearly, for a uniformly chosen permutation tt, E{tt) = 1/n. On the other 
hand, P(c7t) > P(Geom(l/n) > t) = {\ - 1/n)* and thus 

liminf max dc, cAn^ogk + an) > e~^°^^~°^ > e""' / K. 



¥{Xt >C)>1- 

> 1 - 

> 1 - 

> 1 - 
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It follows that d{a) defined as 

d(a) := lim max dr, ci (nlosk + an), 

lies somewhere in the shaded region of Figure [2j □ 

d{a) 
1 - 



l/K 



a 
Figure 2: No cutoff for finite k 

We consider now the random-to-random transposition shuffle. We obtain cutoff by improving 
the upper bound on the mixing time of 1 card and again using coupon-collector arguments 
for the lower bound. 

Proof of Theorem^ We show that the following hold: 

(i) . limsup max dd,. . .,0^ (0.5n log + an) < e~^°. 

n— !>oo (ci,...,Cfc)ef2j. 

(ii) . Suppose A; — > cxo as n ^ oo. Then 

lim liminf max dc, ci, (0.5n log A; + an) = 1. 

a->—oc n^oo (ci,...,Cfe)enfc 

(iii) . Suppose k is bounded above by constant K for all n. Then 

liminf max dc, ci. (0.5n log A; + an) > e~^"/i^. 

For (i) by using similar arguments to the proof in the last lemma, it suffices to show that 
maxjg[„] dj(t) < e~^*(^~^/'^)/". To do this we adapt the proof of Lemma [H we couple now 
both the choices of the left hand and the right hand so that the locations of card i in the 
two decks will become equal the first time that either Lt or Rt chooses card i in iit and the 
other chooses a location other than vr^^ or o"^^. Note that if both hands make choices from 
{a^_} (i) , TTfTj: (i)} the cards will not become matched. We thus have 

P(cj-_i(i) = Tr^_}{i)) > 1 - exp(-2t(l - 2/n)/n). 

For (ii) we use the same argument as in the proof of (ii) in the previous lemma, but with 
a few modifications. We now set T = mi{t > : \{TZt U Ct) H C| = k — C}, where Ct = 
{cri_(Li), . . . , at-{Lt)} is the set of labels of cards chosen by the left hand by time t. We now 



e-"/K 
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obtain E{Xt) = k{l - l/n)^* and E{Xf) = k{l - l/n)^* + k{k - 1)(1 - 2/n)2*. We thus have 
¥{Xt > C) > 1 — e/2 for ah n and then a sufficiently large, with t = O.bnlogk — an. This 
complete the proof of cutoff of the partial mixing time of the random-to-random transposition 
shuffle. 

For statement (iii), we again let E{fj,) be the event that special card ci is in location cJq"^(ci) 
in permutation /x. As before we have -E'(vr) = 1/n. On the other hand, 

P(cTt) > P(Geom(l/n) > = (1 - l/nf\ 

and thus 

liminf max dc, ci, (0.5n log A; + an) > e 



4 Cyclic-to-random upper bound 

We now consider the mixing time of one card in a deck which evolves by the cyclic-to-random 
transposition shuffle. Using the coupling technique for the choices of the right hand from 
Lemma [H when card i is chosen in deck (vr^) with the right hand, we select card i in deck 
{at) with the right hand so that after this transposition cards i will be in the same location 
in their respective decks. We shall also make other modifications to the choice of the right 
hand to speed up the time to couple the cards with label i (we shall refer to them as the i 
cards) . 

We denote by Rf and RJ the locations chosen by the right hand at time t in decks (at) 
and {iTt), respectively. The left hand chooses the same locations in each deck at every time. 
Suppose at time t, Lt selects position (i) (i.e. card i). Further suppose that in a few steps 
at time s the left hand will select location 7r^_}{i). A good choice for R^ would therefore be 
Rf since card i in as will likely still be in position R^ . We formalize this idea in the proof of 
the following lemma. 

Lemma 4. For the cyclic-to-random transposition shuffle, there exists a constant c > such 
that for all n sufficiently large we have 

maxdiit) < cexp(-t/n)((0.237)L0-693t/nJ + l/n). 

Proof We shall use the notation in which a{i) = j says that card j is in position i, so that 
£7(1) is the label of the top card in the deck. Let Ct be the event that the locations of the 
cards with label i in each deck are equal at time t, that is: 

Ct:= {ar\r)=7rr\i)}. 

Let TZt be the event that by time t the right hand has selected card i in vr, that is: 

nt:={3s<t: R:_=7T;}{i)}. 



□ 
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If = Trg}{i) we set Rg = o-g}{i) so that TZt C Ct and therefore 

F{C^) = F{C^\n^t)¥in^t) = F{C^\ exp(-Vn). (8) 

We are interested in matching the locations of the two i cards as quickly as possible. One 
property of the cyclic-to-random shuffle which we shall exploit is that if a card is not selected 
during a round (ra steps) with the right hand, then it will be selected with the left hand at 
least once in that round. 

For each t, we define a random variable 5t which takes value either or 1 and which we 
shall refer to as the phase of the system at time t. At time we set 8q = 1. Let Mt{i) = 
max((j^^(i), 7r^"^(i)) and mt{i) = inin{a^^ {i),TT^^ (i)). We shall define a notion of distance 
denoted Dt between the location of card i in at and the location of card i in vrf. For times t 
of phase 1, we define this distance Dt to be 



Dt := 



n-Mtii) + mtii), if mtii) < Lt < Mt(J,) , 
Mt{i) — mt{i), otherwise. 



Let £ G (0, 1/2) be a constant, which will be chosen later. If at time t, < Dt < en we shall 
say the i cards are close, otherwise they are far. If we are in phase 1 at time t, we will enter 
phase at time 

min{s > t : Lg E {cr("l^(z),7r^^(z)}and0 < Dg < en}. 

This is the first time after time s that the left hand selects an i card when the two i cards 
are close. 

If at time t we are in phase 0, which we entered at time s < t, we will leave it (and return to 
phase 1) at time 

min{r > t : W = fcHi), 7r;i(i)} \ {Lg}}. 
The distance Dt will remain constant during times of phase 0. 

The proof of this lemma uses three different coupling algorithms. The first is a coupling of 
the evolution of the two decks of cards, the second is a coupling of these two decks with a 
certain Markov chain, and the final coupling is of this Markov chain with another, simpler, 
Markov chain. 

We first describe the coupling of the two decks of cards. For each t G N, let Vt be an 
independent Bernoulli random variable with P(Vt = 1) = 1/n. We denote by Tj the time of 
the start of the z*^ phase 0. We shall define a permutation fi^ as 



Tm—1 



T^Tm-^ otherwise. 



We set Urn '■= {'^Tm-iT^Tm) \ l^m- Further, we define 



I m 



and similarly for R^ . The coupling algorithm at time t is as follows: 
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• If there exists s < t such that = 1, choose Rf uniformly on [n] and set = Rf. (In 
this situation we will have already matched the locations of the i cards.) 

• If for all s <t,Vs = 0, but Vt = I, set R^ = tt;'_}{{) and R^ = cj^I_^(i). (This matches 
the locations of the i cards.) 

• If for all s<t,Vs=0, 

— if i ^ {7rt-{Lt),(Jt-{Lt)}, choose independently RJ uniformly on [n] \ vrjl^(i) and 
R" uniformly on [n] \ a^J:{i). 

— if i G {7Tt-{Lt),crt_{Lt)} and Dt = 0, set R^ = RJ to be chosen uniformly on 

— if i € {iTt- {Lt),at- (Lt)}, Dt ^ and there does not exist an m with t = Tm + , 
we choose independently R^ uniformly on [n] \7rt_}{i) and R" uniformly on [n] \ 
a^_}{i). However, if in fact there exists rh with t = Tfn and R^ = v^{i), we shall 
say that the coupling fails and terminate. 

— if i E {vTt- {Lt), at- {Lt)}, Dt ^ and there exists m with t = Tm + D^^ , we choose 
Rt uniformly on [n] \ Rt^ and then toss a coin which lands heads with probability 
l/(n — 1). If it lands heads we say the coupling fails and terminate. If it lands 
tails we set R^ = Rrm- 

Note that if at any time we have Dt = then Dg = for all s > t. Figure [3] shows a possible 
situation in which the locations of the i cards become matched at the end of a phase 0. 



) 



Figure 3: In this diagram showing how cards can become matched, position 1 is the top of 
the deck and n is the bottom 

It is clear that the distribution of RJ is independent of {Rf, ■ ■ ■ ,Rf_i) and similarly R^ is 
independent of {Rf, ■ ■ ■ , Rt-i)- We now show uniformity. For all times t we have 

nat-{R?) = i)= nn-{Rl) = i)= nvt = l) = Vn. 
For times t ^ {r^ + Dr^,Tm} for any m and j ^ i, 

n<^tAR?) = j) = n^ARt) = j) = {-^^ (^) = i/n. 
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For t 



= Tm for some m, and j ^ {MmHO> i^m^(0} 



n. 



For t 



= Tm + for some m, and j ^ {lJ'm{i),Vm{i)} 



n. 



The reason for having the conditional probabihty in the previous equation is since if there 
exists such a t, we know we have not had either success or failure by this time {Bf^^^ = IJ-'^{i) 
iff = 1 which implies coupling success and Rr^ = i^"^ (i) implies coupling failure) . 

For a, 6 G {1, ... , n}, define 



We note that with this definition \\a — b\\ = n — \\b — a\\. 
Let 

Tcouple = min{min{t > 1 : Vt = l},min{t > 1 : 3m:Tm<t and - L^„|| > -Dr„}}, 

be the first time the locations of the i cards become matched. We now define the coupling of 
the two decks with a 3-state Markov chain, denoted {Xt)t>o- This coupling uses the {R")t>i 
and )t>i which we have just constructed. This discrete-time Markov chain has state 
space {C, -F, S}. This will correspond to the i cards being close, far and coupled (success 
state) respectively. We shall construct the chain so that one step of it corresponds to at most 
(1 + £)n steps of the card shuffling. Our interest is in showing that 



for some absorbing state S. We note that we only require coupling the evolution of the decks 
with a Markov chain for times t with event holding. Therefore in the following discussion 
we assume that this event holds for all times referred to. 

We begin the coupling by constructing a sequence of independent random variables denoted 
{Um)m>ii each uniform on {1, ... , en}. The desired property of these random variables is that 
for each m > 1 with Tm < Tcoupiei Um > Dj-^, almost surely. 

If 3£ < m such that WRri — \\ > Dt-^, we simply choose Um uniformly on {1, . . . , en}. 
Conditionally on Tm < Tcoupie; we define three disjoint events, denoted Ej^^, E'^, and E^: 



• -^m := 3 s G [Tm + - \\,rm + WRrn, - II + (1 " s)^), such that IJ,s{Ls) = i, 
Rs G [R^^ — en, Rrm)- We denote by Sm this value of s. 

• E"^ := 3 s e [Tm + WRrm - h ^m + \\Rt^ - L^^ || + (1 - £)n), such that Hs{Ls) = i, 
R^ G [i^^t^,min(-RL + s - Tm - ||-RL - ^r^hR^m + We denote by Sm this value 
of s. 




a ~ b, if a > 6; 
n — b + a, otherwise. 



P(Ct=|7^,C)<P(XLt/(l+,)„J T^S) 



17 



Figure 4: In these diagrams showing card trajectories, the left hand cycles clockwise around 
the circles 



Figure m shows possible trajectories for events and (for event E^ we show two possible 
trajectories for one of the i cards). 

If, conditionally on Tm < ^couple) event E^ holds for some Sm, set Um+i = WRt,^ ~ II = 
Dt^^i- If however, event E'^ holds we set Um+i equal to \\Rsm ~R'rm II with probability a/ (en) 
and with probability 1 — a /{en) we choose Um+i uniformly from the set {a+l, . . . , en}, where 
a = min(||Ls„ - iZ^^^ ||, en). If E'^ holds we set Um+i = Dr^^^. 

We have to check that under this construction of the C/m, they are all independent, uniform 
on [en], and satisfy Um > i?T,„ almost surely. For the independence, we note that the con- 
struction of Um+i is independent of any of the previous distances, Di, . . . , Dm, and therefore 
is independent of C/i, . . . , Um- If events Em or Em hold, it is clear that Um+i is uniform on 
[en] (for event Em, this follows from Rs^ being uniform on [Rr^ — en,Rr^)). If event Em 
holds, uniformity follows from the fact that {{Rs^ — RrmW is uniform on [a] and symmetry. 
The almost sure inequality Um+i > -^t,„+i is clear since -Dr^+i = ll-Rs^ — -Rr™,ll- 

Next we note that the probability of event E^ U E'^ does not depend on the value of . 
We calculate the probability of this union. 

Let ps = ^{Em U Em \ns{Ls) = i). We are interested in calculating y. Note 

that our conditioning on , ,-. , s means that when an i card is selected with the left hand 



it is not able to stay in that position but jumps uniformly to all other n — 1 possible positions. 
For simplicity of notation during this calculation we shall subtract Tm + WRt^ — -^t™,|| from 
the time index. Thus we are interested in pQ. For s G {en,n — en), 



'Tm + (l+en) 




r=s+l 



We solve this using p., 



,_i = 2en/(n — 1), to deduce that for s G (en, n — en). 



Ps = 2e(l - 1/n) 



(9) 
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For s G [0, en], 

^ en n—en—l 

m — 1 n — 1 ^-^ n — 1 ^-^ 

r=s+l r=sn 

Using equation ^ we deduce that for s € [0, en], 

Ps = l + 2e(l - l/n)"+^"-" - (1 - l/n)"-i-=", 

and so, in particular, 

Po = l + 2e(l - 1/n)" - (1 - l/n)-^"-i. 

We now describe how the Markov chain jumps. We shah embed the times of the jumps of 
the Markov chain into the time of the shuffling process. This will enable us to prove (by 
induction) that at time t of the card shuffling, at least [t/{(l + e)n)J jumps of the Markov 
chain have been made. Furthermore, our construction of the coupling (of the Markov chain 
with the card shuffling) ensures that if the locations of the two i cards have not become 
matched by time t then the Markov chain will not be at state S after [t/((l + e)n)J jumps. 
Our embedding of the jump-times is such that the chain can only jump at times t such that 
Lt € {(T^J: (i) , TT^J: (i)} ■ To begin the inductive process, consider the first time the left hand 
selects an i card. If the cards are close at this time, we start the Markov chain from state C 
at this time, otherwise we start the Markov chain from state F. 

Let 

Ml := inf{m > 1 : - || > D^^}, := inf{m > 1 : - || > U„,}. 

We note that if ||i?rm — Lt-^\\ = D^^, we terminate the couplings at time Tm and say they 
have failed. Also, note that Mi < M2 almost surely. Suppose that the chain has just jumped 
to state C (either from state C or state F). This corresponds to a time Tm (start of the m*^ 
phase 0). If we have m = M2, the chain will next jump to state S at time Tm + en. At this 
time the locations of the two i cards will certainly be matched. If m = Mi < M2, we run the 
Markov chain independently from the evolution of the cards from time Tm, according to the 
transition matrix P given below with jumps every (1 + e)n steps and with the first jump (at 
time + (1 + e)n) conditioned on not going to state S. If m < Mi, we are in one of the 
following three situations: 

1. E}^ occurs. Then the next phase starts at time 

Tm+l =Tm+ " \\ + " || < T„ + en + n. 

We make the chain jump to state C at this time. 

2. E"^ occurs. Then the next phase starts at time 

Tm+l =Tm + n+ \\Rf^^ - L^^ \\ < Tm + U + EH. 

We make the chain jump to state C at this time. 
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3. Ef^ occurs. Then the next jump of the chain is at time 

inf{t >Tm + Dr^ : Lt e {at^(i),TT^}(i)}}, 
and it jumps to state F. Note that this t satisfies t < Tm + n. 

In each of these situations the amount of time we have to wait after time Tm until the next 
jump of the Markov chain is less than (1 + e)n. 

On the other hand, suppose the chain has just moved to state F (either from state F or from 
state C). At this time the left hand selects an i card. If it is moved to a location so that the 
two i cards arc close, we wait until the next time an i card is selected (this will be the start 
of the next phase and will be less than n steps later) and at this time make the Markov 
chain jump to state C. On the other hand, if it is not moved to such a location, we make the 
chain jump to state F the next time an i card is selected with the left hand (again this will 
be less than n steps later). This completes the inductive step, showing that each step of the 
Markov chain takes at most (1 + e)n time steps of the shuffling. 

The fact that this does indeed give a process which is Markov follows from the Um being inde- 
pendent, identically distributed random variables and the event E^^ U E"^ being independent 
from the process up to time Tm- 

Consider time r^: the start of the m*^ phase 0, for m < Mi. Using the uniformity of ?7^, 
the probability the chain next jumps to state S is 



en ^ 



en + 1 



x=l 



en ^" ^ 2(n-l)' 



Using po as calculated above we further deduce that (if currently at state C) the probability 
the chain next jumps to state C is 

(1 - (1 - 1/n)— 1 + 2e{l - 1/nr--) . 

This leaves us with a probability of 

((1 - 1/n)— 1 - 2s{l - 1/n)") 

that the chain jumps to F if at state C. 

Let g{n) = 1 — 1/n. We obtain the following transition matrix of the Markov chain {Xm)m>i 
(here the first row/column corresponds to state C, the second to state F and the third to 
state S): 

^(l-5(n)-— i + 2£9(n)") -±i^ (^(n)— ^ - 2e5(n)--") 

2en 1 _ 2£n q 

n—1 n—1 

1 

The probability of the chain jumping from C to S includes the possibility of termination of 
the couplings. We now fix a ^ > and note that for all n > no for some sufficiently large no 
we have both 

+ ' <l + e (10) 



2(n - 1) 2 
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and 



(1 - l/n)-="-i - 2e{l - < (e" - 2ee-'+^) + C (H) 

We define P to be the following matrix: 

/ f (1 - + 2ee^-' + f (e" - 2eei-") + ^(1 - e/2) 1 - f - C 
P ■■= \ 2e 1 - 2e 

\ 1 

Note that we have P{F,F) < P{F,F) and P{F,C) > P(F,C). Furthermore, conditions 
(UnD and (HI]) imply that P{C,C) > P{C,C), P{C,S) > P{C,S) and P{C,F) < P{C,F). 
Therefore, performing the obvious coupling between Markov chain {Xm)m>i and a Markov 
chain {Xm)m>i with transition matrix P we have for all n > no, and each m, ¥{Xm S) < 
P{Xm S). At this point we choose the value of e which minimizes the second largest 
eigenvalue of matrix P (we take ^ arbitrarily close to 0). We find (numerically) that this 
optimal value is about 0.442. With this value of e we obtain the value A of the second largest 
eigenvalue to be approximately 0.237. We deduce that there exists a constant k such that for 
all n sufficiently large (and regardless of starting location), ¥{Xm S) < k(0.237)™. Finally 
we calculate the probability of termination occurring before we couple the locations of the 
two i cards. Note that termination occurs with probability 2/n for every phase that we 
encounter. On the other hand, we couple the locations of the i cards with probability 1 — e 
during each phase 0. Therefore the probability that termination occurs before the couplings 
are successful is 2/(n(l — e)). 

We therefore obtain that F{C^\n^) < c((0.237)L*/i-442"J + i/^) for some constant c and thus 
using equation ([8]) we have P(Cp) < ce"*/"((0.237) L0-693*/«J + i/^), as required. □ 

Remark 2. There are several places in which the coupling argument could be strengthened 
(although this would result in the dynamics becoming more complicated). One of these ways 
would be to improve on the amount of time we wait between successive times of going from 
state F to state F. Indeed, it is clear that n steps is all that is required rather than (1 + e)n 
(since within n steps the left hand will select an i card). 

Proof of Theorem O We prove this result by combining Lemma H] and Theorem [TJ Taking 
6 = 1/4 - e~^/^, we have t^i^{l/4:) < tmix(^~^^^/^)- We wish to find the smallest t such that 
ce-*/''((0.237)L0-693V"J + i/n) < e'^/'^/k. Solving this, we find t « 0.5005n(log /c + C), with 
C = — log(e~^/^/c — 1), which completes the proof. □ 

5 Further work and open problems 

There are several natural further questions we can ask in relation to partial mixing. The most 
interesting question is the following: Does there exist a semi-random transposition shuffle with 
fc-partial mixing time of an log A; (for some A;<n, fc— >ooasn— )• oo) but with mixing time 
/3n log n for some /3 7^ a? If this is not the case, then is there a general way to extend partial 
mixing to the full mixing for any semi-random transposition shuffle? Indeed, it is possible 
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to provide b ounds on the par tial mixing time for larger values of k using the card marking 
techniques of Matthews 1988l | . 



With regards to lower bounds, it can be shown that the partial mixing time of the cyclic- 
to-random transposition shuffle is 0( nlogfc) by obtainin g a lower bound of approximately 



0.12?7-logA; by adapting a method of iMossel et alj 200J]. Determining if and when cutoff 



occurs for both partial and full mixing remains a challenge for this shuffle. 

The method we have developed can be applied to other processes not considered here. These 
include shuffling by fc- cycles, which is known to have a mixing time of (n/k) logn (shown by 
Berestvcki et all |201ll ]) as well as general interchange processes on graphs. Indeed the top-to- 
random transposition shuffle is simply the interchange process on the star graph on n vertices, 
and the random-to-random transposition shuffle is the interchange process on the complete 
graph on n vertices. Studying the fc-partial mixing time is then equivalent to placing just k 
walkers onto different vertices of the graph and asking how long until they reach equilibrium. 
There exists a class of graphs such that for k = o((n/ log n)^/^) we only need to calculate the 
mixing time of one walker on this graph to determine the /c-partial mixing time (with the star 
and complete graphs being in this class). 

Acknowledgements: The author wishes to thank Nathanael Berestycki, Yuval Peres and Perla 
Sousi for useful discussions. 
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